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Abstract 

The one-body density matrix is derived within the Extended Thomas-Fermi approxi- 
mation. This has been done starting from the Wigner-Kirkwood distribution frmction for 
a non-local single-particle potential. The links between this new approach to the density 
matrix with former ones available in the literature are widely discussed. The semiclassi- 
cal Hartree-Fock energy at Extended Thomas-Fermi level is also obtained in the case of a 
S . non-local one-body Hamiltonian. Numerical applications are performed using the Gogny 

and Brink-Boeker effective interactions. The semiclassical binding energies and root mean 
square radii are compared with the fully quantal ones and with those obtained using the 
Strutinsky averaged method. 



1 Introduction 

The one-body density matrix (DM) /j(r, r') = J2a 4'*ai''")4'a{f') or equivalently its Wigner trans- 
form the distribution function f{R,p) (defined below), plays a crucial role in the Hartree-Fock 
(HF) calculations. If zero-range Skyrme forces M are used, only the diagonal part of the DM is 
needed. However, full knowledge of p{r, r') (or f{R,p)) is necessary if one considers finite-range 
effective nuclear forces which are derived from G-matrix calculations in nuclear matter through 
the local density approximation |2|, ^, |^, ^ or postulated empirically with their parameters fitted 
to reproduce some properties of nuclear matter and finite nuclei |^, . 

The full calculation of the density matrix (or the distribution function) is not an easy task 
and requires some computational effort |Q, q|. Consequently, approximations which simplify the 
calculation and, at the same time, show more clearly the physical content of the DM are in order. 
The simplest one is to replace the non-diagonal part of the DM by its value in nuclear matter 
(Slater approach). Finite size effects are added using the density matrix expansion (DME), 
either that due to Negele and Vautherin (NV) |^, [l^ or the modified expansion due to Campi 
and Bouyssy (CB) [^]. Very recently, the CB approach has been applied to HF calculations of 
finite nuclei |1^ using a density-dependent version of the M3Y interaction Q, |5| . 



On the other hand, semiclassical methods |13] are very useful for describing nuclear properties 



of a global character such as binding energies or nuclear densities and their moments. Concerning 
the nuclear ground state properties at HF level, semiclassical approaches are based on the 
Wigner-Kirkwood (WK) ?l-expansion of the distribution function which for a set of nucleons 
moving in a local external potential V{r) up to second order is given by p!^ : 

fwK{R,p) = Q{X-Hw)-—AV6'{X-Hw) 

om 

+ ^[{VVf + ^{p.VfV]6"{X-Hw) + 0{h'), (1) 



where A is the chemical potential and Hy\r is the Wigner transform |13] of the one-body Hamil- 
tonian, which reads 

H^ = ^ + ViR). (2) 

The semiclassical distribution function f{R,p) is a representation of the true phase-space func- 
tion in terms of distributions and is very efficient in order to obtain semiclassical expectation 
values by integrals over the whole phase-space O, |l5|. 



The main purpose of this paper is to derive the exphcit expression of the DM in the Ex- 



tended Thomas-Fermi (ETF) approximation |16] starting from the very recently presented WK 



expansion up to fi order of the distribution function for non-local potentials ||T5|. On one hand, 
we want to establish a link between the NV (and CB) expansions of the DM with this semi- 
classical approach and, on the other, to apply this ETF DM to derive the exchange HF energy 
when finite range forces are used. The paper is organized as follows: In the first section we 
compare the semiclassical ETF density matrix with the former approximations of NV and CB 
in the case of a local potential. In the second section we derive the density matrix and the 
HF energy in the ETF approximation for a non-local potential. We also perform restricted HF 
variational calculations for some selected spherical nuclei using the Gogny Q and Brink-Boeker 
|n] effective forces. We compare these HF ETF results with those obtained quantally, with those 



obtained with the Strutinsky average method |19] and with those which result from the NV and 



CB approaches to the DM. In the last section we give our conclusions and outlook. Technical 
details concerning the calculation of the DM and HF energy in the ETF approach for a non-local 
potential are given in the Appendix. 

2 Extended Thomas-Fermi Density Matrix 

The first step is to perform the inverse Wigner transform of (||) to obtain the semiclassical WK 
density matrix in coordinate space. The definition used here for the Wigner transform of the 
one-body density matrix is [13|: 



f{R,p) = jdse-^P'l^p{R+'-,R-'-\ (3) 

where R = {r\ + r'2)/2, s = r\ — r2 and p are, respectively, the centre-of-mass, the relative 
coordinates and the phase-space momentum. 

After some lengthy but straightforward algebra the semiclassical DM in terms of R and s 
at WK level is given by: 



+ T^ — 7 [joikps) - Akpsjiikps) + kps'^J2{kFs)] 



-rl^^^[kF^kF-]-[ - ^kpshikps) + kls^J2{kFs)\ + 0{n^), (4) 

487r^ kp s s 



where kp = \/^{X — V{R.)) is the local Fermi momentum, ji{kps) are the spherical Bessel 
functions and g stands for the degeneracy. 

This expression, although written in a slightly different way, coincides with the ones obtained 



previously by Dreizler and Gross [0] and Jennings |18|. The first term of the expansion (| 



corresponds to the Slater approach, whereas the h terms are the part that take into account 
quantal finite-size effects. 

The WK density matrix in coordinate space depends on the angle between R and s, however 
for practical purposes and following previous literature [^, 11, 10 1 we perform the angular average 
of eq. (^ obtaining: 

p{R,s) = pq{R, s) + p2{R, s) 

_ gk%?,ji{kps) g ., r„.,, . „, . ., ^ ,22 



+ -—-^/^kp[QJQ{kps) - 3kpsji{kps) - kps J2{kps)] 
[3jo{kps) - 9kpsji{kps) + 2k%s^J2{kps)]. (5) 



67r2 kps 1447r2 



1447r2 kp 

The diagonal part (s=0) of eq.(|5|) is the well-known WK expression for the local density (with 
degeneracy g) |[l^, |l5| : 

To obtain the DM in the ETF approach we shall express the Fermi momentum and its 
derivatives in terms of the local density and its gradients. First, the local Fermi momentum is 
obtained by inverting eq.(y): 

^F = fco-^[^3- + 2^]> (7) 

where k^ = {67r'^p/g)^'^. Notice that for inverting the gradient terms in eq.(^ it is enough to 
replace kp by ko to be consistent with the 7i-order in the expansion of the Fermi momentum 
(R). Writing the gradients of /cq in terms of the spatial derivatives of the local density 

.„,Vt,)^ = f!(^ (8) 

35 p 
g -J p 



one finally obtains the Fermi momentum as: 

,, = (5^)1/3 + i.(5!i!^)-./3[(Zrf!_ 2^], (10) 

g 72 g p^ p 

where the first term of the right-hand side is the pure Thomas-Fermi part and the second term, 
which contains derivatives of the local density, is the % contribution. The semiclassical density 
matrix for a local potential in the ETF approach is obtained from (||) by expanding consistently 
the Fermi momentum kp up to % -order with the help of eqs.(^ - [lO| ): 

Petf{R, s) = p — + —Apijoikps) - 6— 

kps 72 kps 

'' ^^''^%Mk..)-9il^], (11) 



216 p kps 

where now kp = ko = {67r'^p/g)^'^. 

Let us now analyze the main properties of this semiclassical approach as compared with 
the quantal case. Following refs.[^, y, the quantal DM averaged over the s direction can be 
approximated by: 

p{R,s) = l-Jdne''P/^p{r,rX=r'=R = Msl)p{r,rX=r'=R^ (12) 

where p = —ifiiWi — V2)/2 is the relative momentum operator. Expanding the Bessel function 
in a Taylor series one gets: 

n=0 ^ ' n=0 ^ ' 

where M2n are the momentum weighted integrals defined as ||8|: 

^2n = J^^jMlrf{R,P) = irPiry)\r=r'=R- (14) 

It should be pointed out that series (13) as it stands is not useful because it converges very 
slowly and cannot be truncated for large s-values if the even moments are different from zero. 

At this point there are two possibilities for approximating the exact DM. One is to sum the 
series (|l^ ) using some approach for evaluating the momentum weighted integrals and the other 
is to rearrange the series (13) in such a way that truncation is possible. 



First of all, we will show that the semiclassical ETF approach to the DM eq.( O) corresponds 



to the whole sum of the series eq.(13) if the moments M2n are calculated in the ETF approxima- 



tion. Starting from the semiclassical distribution function, eq.(|T|), the WK momentum weighted 



integrals up to fi -order are easily derived. Expanding consistently the Fermi momentum with 



the help of (10) and using eqs.(|8[^ one finally obtains the ETF weighted integrals in terms of 



the local density and its gradients: 

where again kp = ko = {Qir'^p/gY'^ . 

The first term of the right-hand side of (p^ is the Thomas-Fermi weighted integral while the 



second term is just the Ti correction in the ETF approach. For n = and n = 2 one obtains: 

M^^^ = p (16) 

M-.|p./ + l(Zf:l! + lA„ (17) 

which are just the semiclassical counterparts (at ETF-?i level) of the zeroth and second-order 
quantal momentum weighted integrals: Mq = p and M2 = r — A/9/4, where r is the kinetic 
energy density. Notice that in the ETF-?i approach only second-order gradients of the local 
density appear in M^j^^ for any value of n. However, higher order derivatives will appear in 
the moments if the ETF expansion is pushed to higher powers in h. 

Taking into account the Taylor expansion of the Bessel functions in eq. ( [Tl| ) and after some 
algebra one finds: 

rr. ^ ^jlikps) ^ (-If s2"+2 2nr8^ + 3(V/j)2 2n - 3 ^ , 

^--^^'^) = ^4^-SVT3)r(-+^)^^ [ 108 p -^^^^] 



n=0 



2n-F3)(2n + l)! ^^{2n + l)\' ^" 2n-h3' 



oo (_i)ng2n 



= y )-'> \, MfJ^. (18) 

^Q(2n + 1)! 2" ^ J 

From this result it is clear that in the ETF approximation to the DM all the momentum weighted 



integrals appearing in eq.(12), evaluated within the same semiclassical approach, are consistently 
summed. 



Another possibility for approximating the quantal DM is to rearrange the terms in eq.([l^) 
in such a way that the leading term is the Slater term. This is, actually, the way in which the 
NV and CB approaches to the DM are done. For the sake of completeness we shall once again 
briefly derive the NV and CB approaches to the DM following the method outlined in ref. H. 



The starting point is the the identity 

Mab) = £(-ir(4n + 3)^?I^^^^±i^ (19) 

n=0 "''^ °/^ 

valid for any k such that — 1 < b/k < 1 and where P2n+i{x) are the Legendre polynomials 

P2n+l{x) = E™=oCa^'^"""^+' With 

(-1)™ (4n-2m + 2)! 



22n+i rnl{2n - m + l)!(2n - 2m + 1)! ' 



(20) 



Using eq.(|l9|) in eg. ([12D , the angular averaged DM eq.(|l3| ) is also written as: 



^^ 2n + l ! ^„ aj? ^^'^ rn) , k J 

n=0 ^ ' m=0 u 

where j^n+iC^s) = [An + 3)\lJ2n+i{ks) / (ks)'^"''^^ are the spherical Bessel functions normalized to 
unity at s = 0. 

Of course, the semiclassical petf{R, s) (|18|) also fulfills eg. (pi]) . Starting from this eguation 
if the moments are obtained in the ETF approximation ([l5|), the Bessel functions J2n+i{ks) are 
expanded in a Taylor series and the even powers of s are properly sorted out, one recovers (|lq). 

The NV approach consists of keeping the two first terms of the expansion (^) and taking 
k = kp. In this case the DM reads: 

/„ N ^iiikps) 35 i-jf/cirs) .3 , 9 . ^ ,^, 



The CB approximation keeps only the first term of (|2l| ) but with k fixed in such a way that the 
second term of (pTI) identically vanishes: 

PCB[R, s) = p — ~ , (23) 

kps 



where kp = \/5M2/3p. 

The NV and CB approaches are truncations of the true expansion of the guantal DM eg. ([l^) . 
As is discussed in |^ , in these approximations only the Mq and M2 momentum weighted integrals 
correspond to those obtained with the exact DM, whereas any higher momentum weighted 
integral in these approaches M2X (A > 1) has little to do with its exact guantal value. 



Let us now discuss the results obtained using the different approaches to the DM analyzed 
previously. To do this we consider a ^^Ca nucleus described using harmonic oscillator (HO) 



wavefunctions with an oscillator parameter a = ^/rnujjh = 0.516 fm~^. Figure 1 displays the 
ratio of the off-diagonal to diagonal DM p{R, s)/p{R) as a function of the interparticle distance 
s for selected values of the distance from the centre-of-mass R. The different curves shown in 
this figure correspond to the quantal DM (black dots) , the NV (dashed-dotted line) and CB 
(dashed line) approximations and the semiclassical ETF approach calculated using the quantal 
local density (solid line). As general trends the approximations to the DM considered here 
reproduce reasonably well the quantal values in the range of i? = 1 — 3 fm, they show some 
deficiencies at i? = fm and clearly start to separate from the quantal values for i? > 4 fm. 
However, in the whole range of -R- values analyzed, the quantal approaches NV and CB reproduce 
the quantal DM better than the semiclassical ETF calculation. 

At this point two comments are in order. First of all, it should be pointed out that all the 
approaches to the DM considered in this paper are, in fact, distributions (see eq.(||) for the 
ETF DM and ref. ||8| for the discussion of the NV and CB cases) and consequently, the only 
meaningful comparison should be done through the moments in k and R spaces. 

On the other hand, the semiclassical approaches to the DM are obtained by switching off 
shell effects. Consequently, the Slater and ETF approximations to the DM should be compared 



with the smoothed DM obtained using the Strutinsky averaged occupation numbers |19| rather 
than with the quantal DM. To do this, one starts from the smooth distribution function, which 
for closed HO shells reads j20|: 

oo 

fstiR,p) = 8j2i-lfhKe-'Lli2e), (24) 

K 

where e = p^ /m + muj^B?', L^ are the generalized Laguerre polynomials and fix the Strutin- 
sky occupation numbers [^]. Performing the inverse Wigner transform of (|2J) and averaging 
over the angles, one obtains the Strutinsky DM in coordinate space to which the semiclassical 
approximations (Slater and ETF) should be compared. 



^^*(^'^) = 7(2^^^*^^'^)^^'''''' 



3 oo o J^ 2 2 



E nKe-^(^^-x) E i-l)-^Lf^i2a-R-)Lf_,s'^). (25) 



7r3/2 

K=0 Kx=0 



The NV and CB expansions of the DM can also be considered within the semiclassical framework 



if the M2 moment is calculated in the ETF approach. From (18), it is clear that in this case NV 
and CB become truncations of the ETF density matrix. It is interesting to look at the quality 
of these approximations because they have been used for calculating the exchange part of the 



nucleus-nucleus potential (see |21| and references quoted therein). 

The Strutinsky DM (|25| ) for ^^Ca is obtained with a HO parameter a=0.516 fm^^ and with 
a smoothing parameter [|l^, |l^ 7=1.25 hto. Figure 2 collects the semiclassical results where 
the ratio p{R,s)/p{R) is displayed for the Strutinsky smoothed DM (|25| ) (black dots) and for 
the ETF DM (solid line) as well as for the semiclassical NV (dashed line) and CB (dashed- 
dotted line) truncations of the ETF DM. Notice that in order to completely remove the shell 
effects, the ETF, NV and CB density matrices have to be obtained using the Strutinsky local 
density [^]. From this Figure it can be seen that the ETF ratio reproduces reasonably well 
the Strutinsky result. The ETF quotient is better than the NV result in the whole range of R 
and s distances analyzed and better than the CB for small values of R. The difference between 
ETF and Strutinsky ratios indicates that the h- expansion in ETF does not fully converged. 
Consequently, it would be necessary to add h contributions to the ETF DM to obtain ETF 
expectation values in good agreement with the corresponding Strutinsky results. 

3 Restricted variational energy calculations 

The second part of this paper is devoted to discussion of the ability of the ETF approach to 
the DM for describing the HE binding energy of finite nuclei when effective finite-range nuclear 
interactions are considered. 

First, we derive the ETF approximation to the DM starting from the non-local HE potential: 

V"P{r, r') = V^{r, r')S{r - r') - V^{r, r'), (26) 

where V^ and V^ are the direct and exchange parts of the HE potential. In Wigner represen- 
tation eq.(P6|) becomes: 

V{R, k) = VdiriR) - VeAR, k), (27) 

where 

VdiriR) = J dR'v{R, R')p{R') (28) 



and 

VeAR,k)=gj j^wik,k')f{R,k'). (29) 

In these equations v{R, R') is the two-body effective interaction, w{k, k') is its Fourier transform 
and g stands for the degeneracy (for the sake of siniphcity we consider a simple Wigner force in 



eqs.(28) and ([29|)). Consequently, the Wigner transform of the one-body Hamiltonian will be; 



Hw{R,k) = %^ + V{R,k). (30) 

If the HF potential is spherically symmetric in k, i.e. V{R,k), the WK distribution function 
required for semiclassical calculations is |15|; 

fwK{R,k) = — fwK{R,k)dnk 

= e{X-Hw)-^S'iX-Hw)FiiR,k) + ^6"{X-Hw)F2{R,k), (31) 

0771 24:771 

where the functions Fi{R, k) and F2{R, k) are given by: 

F,iR, k) = ^[^AViSf + kfk) - k\Vff] (32) 

F2{R, k) = _[-^(Vy)2(3/ + kfk) + k^fAV - 2fc2/Vy V/]. (33) 

In eqs.(32- 33), / is the inverse of the position and momentum dependent effective mass: 

777 77? 

m*{R,k) h^k 

where the suscript k indicates a partial derivative with repect to k. 

Due to the fact that the effective-mass corrections are included in the h part of the distri- 
bution function (|3l]), they are calculated using the h order exchange potential in eq.(p^) to be 
consistent with the 7i-order in the expansion of the WK distribution function. 

Following the steps indicated in the Appendix, the ETF density matrix for each kind of 
nucleon in the case of a non-local potential can be written as: 

- [(18 + QkA)t^ _ 2,Jo{kFs)\Ap 
f kps 



/ f J J f^ kps 



10 



where now kp = (Svr^p)^'^ and the inverse effective mass / (^) and its derivatives are computed 
at k = kp. We use here g=2 because in this way the ETF DM eq.(p5|) can be directly apphed 
to non-symmetric nuclei. 

If all the spatial and momentum derivatives of the inverse effective mass are dropped in 



eq.(35) one recovers the ETF DM for the local case eq.(ll). If only the momentum derivatives 
of / are eliminated, one obtains the ETF DM corresponding to the case of a position-dependent 
effective mass. This latter case is just the situation that appears when one uses zero-range forces 
such as the Skyrme interactions. 

The next step is to obtain the ETF approach to the HE energy, which for an uncharged and 
spin-saturated nucleus can be written as: 



Ehf = Y.J<^^ ^^^ + lp{R)Vdir{R)-\JdsVe^{R,s)p{R + ^,R 



(36) 



where the subindex q stands for each kind of nucleon. 

The HE energy in the ETF approximation is obtained as explained in the Appendix and 
reads: 

h^ ._. 1 



Eetf = Y.J'^^ 



^^TETFiR) + ^p{R)Vdir{R) - eg^^iR) 



(37) 



In this equation tetf{R) is the kinetic energy up to h order for each kind of nucleon in the 
ETF approximation and reads 

TeTf{R) = TETFfi[R) +TeTF,2{R) 

= -^kEp + ^^^[1 + -^kFj + -kE— - -kEj,] + -Ap[4 + -kEj] 

and e^J^ {R) is the exchange energy density for each kind of nucleon in the same approximation 
given by 

+ |1|(/ - l)(r..P - Iklp - iAp) + *=pA(^^ - 1; Artl, (39) 

where v{s) is the central nucleon-nucleon interaction, kE = (Stt^p)^'^ and / and its k derivatives 
are calculated at k = kE- For Gaussian type forces such as the Gogny or Brink-Boeker interac- 



11 



tions used in this paper, the exphcit form of the lowest order exchange energy £^J'q{R) can be 
found in ref. [p^ j. 

In the special case of a local potential the fi^ part of the kinetic energy density eq. (|3q) reduces 
to the well-known result 

rETF,2 = ~^^ + \^P. (40) 

6b p 6 

if the effective mass is not momentum dependent, one recovers the result of IT 

^(Vp)2 1 IpA/ I VpV/ 1 pjVff 

36 p 3^6/6 / 12 p 



TETF,2 = — : h-A/>+-— ^ + - J — -2 . (41) 



For the particular case of the Coulomb potential, the direct calculation of the exchange energy 
density up to fi^ order in the ETF approach (local case) leads to 

ETF ^_3/3.i/3 4/3 ^ (yp)'^ (A0-] 

^Coul,ex ^{^) P 4327r(37r2)l/3 p4/3 ^^ ) 



where p is the proton density. Eq.(|42D agrees with the result reported previously in [17|. 

As a first test of our ETF approach, let us compare the exchange Coulomb energy obtained 
using eq.(E^) with the quantal result as well as the same energy derived through the NV, CB and 



Slater approximations to the DM. To this end and following ref. [11|, we use HO wavefunctions 
with fixed parameters a = 0.752fm for '^He, 0.546 fm for ^^O and 0.481 fm for '^^Ca. Table 
1 shows the quantal (QM label), NV, CB, Slater (SL label) and ETF results for the exchange 
Coulomb energy. From this Table it can be seen that the ETF results almost reproduce the 
quantal values and improve those obtained using the NV, CB and Slater approximations. 

3.1 Comparison with quantal results 

In order to check the quality of our approach in the calculation of HE energies, we present here 
a restricted variational calculation for uncharged ^He, ^^O and '^^Ca nuclei using the Brink- 
Boeker and Gogny forces. To this end we use HO wavefunctions and minimize with respect to 



the HO parameter a = \Jmu:l'h. 

At this point it should be noted that the semiclassical energy (B^) (as well as that calculated 
in the simplest Slater approach) is free of shell effects if it is obtained using a smooth density [^] . 
Consequently, the semiclassical results should be compared rather with a HE calculation based 
on the smoothed Strutinsky density matrix (p^) than with the purely quantal HE calculations. 

12 



However, a direct comparison with quantal results in the mean field approach is possible if 
shell effects are added to the semiclassical results (at Slater or ETF levels) according to the 
Strutinsky energy theorem. One possible way of incorporating shell-effects is based on the 
Kohn-Sham p4l approximation widely used in atomic physics ||T7[| and discussed for the nuclear 



case in [25|. Basically the KS approach consists of solving the quantal HF equation using a 
local effective potential obtained as a functional derivative of the density-dependent exchange- 
correlation energy. In our restricted variational calculation the Kohn-Sham scheme is applied 
by minimizing the sum of the quantal kinetic plus direct energies with the semiclassical (Slater 
or ETF) exchange energy eqs.(p9|). 

The exact quantal energies (corrected from the centre-of-mass motion) are obtained from 
eq.(|36[) with the DM evaluated analytically ||2^ (however, see below for the special case of closed 
HO shells) together with its corresponding root mean square radius (RMSR) < r^ >^'^. These 
quantities are displayed in Table 2 with the label QM for the Brink-Boeker and Gogny force. 
Table 2 also shows the same results obtained using the NV and CB truncations of the quantal 
DM (labelled NV and CB respectively). The binding energies and RMSR obtained within the 
KS scheme starting from the semiclassical Slater or ETF exchange energies eq.([3^) are also 
collected in Table 2 with the SL(KS) and ETF(KS) labels. The differences between the KS 
(Slater and ETF) results and the purely quantal ones show the quality of the semiclassical 
approach to the exchange energy. From this comparison one can see that the Slater approach 
is very poor in the case of the Brink-Boeker force, underbinding all of the considered nuclei 
and giving RMSR larger than the quantal values. However, the result is more satisfactory for 
the Gogny force. This difference is due to the fact that the non-local effects are larger in the 
Brink-Boeker force than in the Gogny case. The non-local effects are better accounted for in 
the ETF(KS) approximation for which agreement with the quantal HF results is good for both 
effective forces considered in this paper. The ETF results in the KS scheme are similar to those 
obtained using the NV and CB approximations to the quantal DM for both Brink-Boeker and 
Gogny interactions. 

As has been pointed out in Section 2, the NV and CB truncations of the quantal DM become 
truncations of the ETF DM if the M2 momentum weighted integrals are also computed with 
the same ETF approach (0). To check the quality of these approximations to the ETF DM, we 
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have again computed the binding energies and RMSR using the NV and CB truncations of the 
ETF DM within the KS scheme. The corresponding results are also collected in Table 2 with 
the NV(KS) and CB(KS) labels. From the analysis of the KS results it can be seen that the 
agreement of NV(KS) and CB(KS) with ETF(KS) is similar to the one found comparing the 
NV and CB results with the QM values. On the other hand, this agreement improves when the 
non-locality of the effective force is smaller. From Table 2 it can also be seen that the ETF(KS) 
energies and RMSR agree better with the corresponding quantal results than the NV(KS) and 
CB(KS) values. 

Some time ago another different approximation was presented in the literature j^^- ^^ 
this approach rather than starting from first principles, a phenomenological density matrix is 
proposed in which the parameters were determined by imposing the correct local semiclassical 
kinetic energy density and the projector character of the DM in an integrated form. Using this 
parametrized DM and the Gogny interaction, the binding energy of ^^O and '^^Ca (including the 
Coulomb energy) are —128.3MeV and —337.9MeV respectively. These results can be compared 
with the ones obtained in our ETF(KS) approach —122.6MeV and —335.2AdeV and with the 
fully quantal results —126.0MeV and —338AAIeV. The results obtained using the ETF(KS) 
approximation are slightly worse than those obtained with the phenomenological DM discussed 
previously, but clearly improve the CB results reported in [p7| , which are —118.2MeV and 
-326AMeV for ^^O and ^°Ca respectively. 

3.2 Comparison with Strutinsky results 

Let us now discuss the quality of the Slater and ETF energies within the semiclassical framework. 
In this case they have to be compared with the ones obtained using the Strutinsky averaged 



method |1£]. The starting point for a Strutinsky calculation of the energy using trial HO 
wavefunctions is the smooth density matrix ( p5[ ) from which the particle and kinetic energy can 
also be derived. For an effective nuclear interaction with two Gaussian type form factors (as 
in the case of the forces studied in this paper) and HO closed shells, the direct and exchange 
energies can be obtained analytically: 



2 "] iiV / ,-,,2 \ 3/2 ^max J^max KM 

E..r = T.-^(^) E E nKUM E E 4/-x.(o)4f-M.(o) 

.•_i V " \ li / i^—n n,f— n ^, — n fl,r, — n 



14 



r(i^i+Mi + |)/a2xi^i+Mi 



Ki\Mi\ 



It 



F 



Mi,-Ki,-Mi-Ki 



a'' 






(43) 



E, 



E 



32X,,/a2\3/2^_ 



max ^^max 



Kmax—KKmax—K T (' ZV' _L 3\ 

^ '* ^ K=Q Ki=Q Mi=0 



K! 



r(iri + Ml + 



KiIAfi! 



1 



7i 



2\ K^+hh 



F 



Mi,-Ki,-Mi-Ki-^,1-^{1 



a'^-,-2 



(44) 



where 7j = 2//i? + a^ and F are the Gauss hypergeometric functions. In eqs.([4^^) X^^i = 
Wi + hi/2 — hi/2 — mi/4 and X^^i = Wi/A + hi/2 — hi/2 — rrii are the usual combination of the 
direct and exchange parameters of the central effective interaction and /ij is the range of each 
Gaussian form factor. 

The Strutinsky occupation numbers that come into in the energy calculation are obtained 
from a HO spectrum. In this way the smooth energy becomes a function of the HO length a. The 
Strutinsky energy is obtained by minimizing with respect to a to simulate the self-consistency 
P9| with the additional constraint that the minimization procedure is performed in the plateau 
region p3| , ll^]. With a smoothing parameter 7=1.25 huj, the HO parameters that mininimize 
the Strutinsky energies of the '^He, ^^O and '^^Ca nuclei are a=0.647, 0.550 and 0.509 fm^^ 
respectively using the Brink-Boeker force and a=0.643, 0.567 and 0.516 fm~^ in the case of 
the Gogny interaction. The binding energies and RMSR obtained in this way for uncharged 
^He, ^^O and '^^Ca nuclei are collected in Table 3 with the label ST. The semiclassical binding 
energies at Slater and ETF levels are computed starting from the Strutinsky particle density 
obtained previously in order to drop the shell effects completely [ p^ . These results are shown 
in Table 3 labelled SL and ETF(a) respectively. 

The Strutinsky value represents the energy which varies smoothly with the number of nu- 
cleons A. For each nucleus the difference between its quantal value QM reported in Table 2 
and the corresponding ST result given in Table 3 is the so-called shell energy. This is a subtle 
quantity that is not reproduced by SL or ETF approaches up to h order. As has been pointed 
out in previous literature, if the ETF kinetic energy density functional only contains the h and 
^2 contributions, its integral is not able to reproduce the Strutinsky kinetic energy at least in 
the case of a set of nucleons moving in a HO or Woods-Saxon type external potential, p^]. 
However, if the h contributions are included in the functional, the ETF kinetic energies are in 
much better agreement with the Strutinsky values p^, |2M|. In our non-local calculations the 
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differences found between Strutinsky and ETF (up to h order) total energies are roughly similar 
to those found for the kinetic energy in the case of an external HO potential |2^. This fact 
suggests including approximately the h corrections to ETF by adding to tetf which enters 
equations (|38|) and (|3^) the h contribution in the local potential case: 



From this approximated calculation we find that almost all the correction comes from the kinetic 
energy term, whereas the exchange part gives only a minor contribution. The total energy when 
this h correction is included perturbatively is shown in Table 3 labelled ETF(?i ). It can be 
seen that the Strutinsky binding energies are very well reproduced by this approximate ETF(7l ) 
calculation showing again the importance of including h corrections in ETF in order to obtain 
a better description of the shell energies [22, 25]. 



For finite-range forces the non-local effects contribute to the DM (35) through the gradients 
and the derivative with respect to k of the inverse effective mass calculated at k = kp = 
{Qn'^p/g)^'^. To investigate the influence of these non-local corrections to the HF energy, we 
have again obtained this energy using the DM corresponding to a local potential (pTJ). In this 
case, the kinetic energy reduces to that corresponding to the local case and the h exchange 
energy is calculated using ( p9D but with r corresponding to the local case. The HF energies 
calculated in this way are also displayed in Table 2 with the label ETF(b). From these results it 
can be seen that, in fact, the h effective mass corrections to the DM ( p5|) are almost negligible 
for the Gogny force but they become more important for the Brink-Boeker interaction where 
the non-local effects are larger. 

4 Summary and Outlook 

In this paper we have derived, for the first time to our knowledge, the Extended Thomas-Fermi 
approximation to the one-body density matrix up to h -order for a non-local single particle 
Hamiltonian. The h contribution can be written in terms of spherical Bessel functions combined 
with second-order gradients of the local density and the inverse of the effective mass as well as 
momentum derivatives of the latter computed at the Fermi momentum. This density matrix 
includes, as particular cases, results reported previously in the literature [O, O] for the local 
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case. 

We have compared this new semiclassical approximation with former approaches, namely the 



Negele-Vautherin ||9[ and Campi-Bouyssy |11] ones. It is found that as in the case of the quantal 
density matrix, the Extended Thomas-Fermi approximation sums all the momentum weighted 
integrals Q, but with their quantal values replaced by their semiclassical counterparts. In this 
respect the Extended Thomas-Fermi approach differs from the Negele-Vautherin and Campi- 
Bouyssy approximations, which are truncations of the quantal density matrix. It should also 
be pointed out that if within the Negele-Vautherin and Campi-Bouyssy approaches the quantal 
momentum weighted integrals are replaced by their semiclassical counterparts, they become 
truncations of the Extended Thomas-Fermi density matrix. 

We have applied this semiclassical approach for deriving the Hartree-Fock energy of a nucleus 
in the case of effective finite-range interactions. In this case the h -order Extended Thomas- 
Fermi kinetic and exchange energies contain, in addition to the second-order gradients of the 
local density and inverse effective mass, new terms that account for the momentum dependence 
of the effective mass. 

We have checked our Extended Thomas-Fermi approach to the Hartree-Fock energy by per- 
forming restricted variational calculations with the Gogny and Brink-Boeker effective forces of 
the binding energy of '^He, ^^O and '^^Ca using trial harmonic oscillator local densities. At 
a quantal level our Extended Thomas-Fermi approach up to h -order to the exchange energy 
can be used within the Kohn-Sham scheme to obtain a local approximation to the non-local 
quantal energy density. It is found that this Kohn-Sham approach gives binding energies and 
root mean square radii similar to the corresponding results obtained using the Negele-Vautherin 
and Campi-Bouyssy truncations of the quantal density matrix. 

We have also performed restricted variational Strutinsky averaged Hartree-Fock calcula- 
tions with the same finite-range nuclear interactions. These results are compared with the 
pure Thomas-Fermi (Slater) and Extended Thomas-Fermi values obtained starting from the 
Strutinsky kinetic energy and particle densities. The Extended Thomas-Fermi binding energies 
including only h -order contributions are not able to reproduce the Strutinsky results and con- 
sequently cannot be used for obtaining the shell energies. We have approximately estimated 
the h -order contribution to the binding energy and verified that if this correction is taken into 
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account, the Strutinsky values are practically recovered. 

Although the Extended Thomas-Fermi approach applied to a non-local one-body Hamilto- 
nian gives reasonably good results, to improve this semiclassical approximation by adding the 
full h -order contributions seems to be in order. Another way of improving the semiclassical 
results presented here to obtain the smooth part of the energy is by using the Variational Wigner- 



Kirkwood approach [15| which slightly differs from the Extended Thomas- Fermi approximation 
presented in this paper. We reserve these extensions of the semiclassical calculations in the 
non-local case for a future work. 

On the other hand, other effective finite-range forces such as M3Y together with the Double 
Folded Model have been recently applied to compute the real part of the microscopic heavy-ion 
optical potential. In these calculations, the exchange part is usually obtained using the Negele- 
Vautherin or Campi-Bouyssy approaches to the density matrix with a semiclassical kinetic energy 



density [21|. To use the full Extended Thomas-Fermi density matrix to obtain the real part of 
the heavy-ion optical potential is another promising application of the method developed in this 
paper and will be presented in a forthcoming publication 
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Appendix 



The WK density matrix up to h order (assuming degeneracy g) in coordinate space for a 
non-local potential is given by the inverse Wigner transform of ( pl|) 



Pwk{R,s) 



+ 



(2^)3 



^^ }wK{R,k)e^^' 



Gvr^ kps 



dkk'^Jo{ks)Fi{R, k)6{X - Hw) 
dkk^Jo{ks)F2{R, k)6{X - Hw) 



that written in terms of the local Fermi momentum kp with the help of: 

m 5{k — kp) 



(46) 



S{X - Hw) 



f^kpfiR^kp) 



and 



reads 



dx_ 

dkp 



f^k 



m 



^+V,{R,kp) = -^'^^'^^^'^ 



m 



(47) 



(48) 



.„ X gkl^3ji{kps) , . ., ^ ,„. 

Pwk{R,s) = 77-7 — ; ^Jo{kps)p2,WK{R) 



67r^ 



+ 



9 djo{kps) 



kps 
m 



487r2 dkp h^kpfiR, kp) 



m 



d , kF2 



+ 



m d 



Fo 



h^kf{R,k)dk'f{R,k)' n^dk'f{R,k) 



3kFi 



f{R,k) 



+ 



g d'^joikps) 



m 



k=kp 487r2 94 'h^kpf{R,kp)^ 



kFo 



f{R,k) 



k=kf 



(49) 



where P2,WK is the WK h -order contribution to the density in the non-local problem given in 



the Appendix A of |15| and the inverse effective mass f{R, k) is defined as in eq.(|3j). 

In eg . (ps|) the gradients of the non-local potential V{R,k) appearing in Fi{R,k), F2{R,k) 
and their momentum derivatives have to be evaluated at k = kp. To do this one starts from the 
definition of the Fermi energy: 



h^kl 

2m 



+ V{R,k)=X, 



(50) 



where kp is also a function of R. Now taking the gradients of (50), the spatial derivatives of 
the potential are transformed into gradients of the local Fermi momentum through: 

n'^kp 



Vy)fc^ + -f{R, kp)Vkp = 



m 



(51) 
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^2 



{AV)k^ + —[kFf{R,kF)AkF + f{R,kF)iVkFr + 2kFVf{R,kF)VkF 



m 



+ kFfkiR,kF)iVkFf]=0. (52) 

To obtain the ETF DM one proceeds as in the local potential case. First the WK local density 



up to ?l^-order is inverted to obtain kF 



2f_ 



kF = ko- —^P2,WK[kf)\, (53) 



where k^ is the zeroth order local Fermi momentum given by k^ = (Gir'^p/g)^'^ and the gradients 
of V that appear in p2,WK [0] have been replaced by gradients of ko with the help of eqs.(51 - 



ID- 

Finally, expanding kF in the WK DM (^) with the help of ( |53|) one obtains the ETF density 
matrix written as 

gko Sji{kos) 



PETF 



Gvr^ k^s 



fk , n,i2fkk j2fk\i'^ko) , f. , „, /) 



fc- 



(7 + 6ko^ + 2k'o^ - kt^y--^ + (4 + 2A:o^)Afco 
J J J f^o J 



gKps 

1447r2 _ 
+ (6-2A:o-r) 7 \-^ko \-2ko—^ fco — 7 -^ — ^ 



/' / " / V " / 



2^2 



1447r2 



(Vfcf 



\2 



A;o 



Afci 



d hikos) 



(54) 



If now the gradients of /cq are written in terms of gradients of the density using eqs.(P -pi) one 
obtains the ETF DM for a non-local potential written as a functional of the local density which 
is just (|35[). 

The HE energy of an uncharged and spin-saturated nucleus in the ETF approximation can 
be obtained by replacing the quantal integrand in (p6|) by its corresponding ETF approximation. 



The ETF kinetic energy density can be derived from the ETF DM using [13| 

tetf{R) = Qa/j - AApetf{R,s)\s=o (55) 

from where eq.(p8D is easily obtained. 

The direct energy is obtained using the the diagonal part of the DM that in the ETF approach 
reduces simply to the local density p. 
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The exchange potential is given by Vex{R,s) = 'v{s)p{R,s) and, consequently, the ETF 
exchange energy will be 



,ETF 



iR,s) = ^JdsV,''/^iR,s)pETFiR,s) 

= 2 / dsv{s) PETFfiiR^ s) + 2pETF,o{R, s)pETF,2{R, s) 



\p\R) J dsvisfJ0^ + J dsVgY{R, S)PETF,2{R. s). 



(56) 



(kEsV 

The integral over s in the h contribution to the ETF exchange energy can be performed 
analytically taking into account the fact that in Wigner space the h ETF exchange potential 
(Slater) can also be written as 



ViyiR,k) = ldsViyiR,s)e''^^ = ldsvis)poiR,s)joiks) 

if the exchange potential is spherically symmetric in k. 

The k derivatives calculated at k = ko are easily obtained starting from (I 



(57) 



KfioSfc(-^'^o) = / dsv{s)po{R,s)s 



dsv{s)po{R,s)s ^^, ^ = \k=ko 



d{ks) m 

^ d^joiks) _ h\f + kfk-l) 
d{ksf 



m 



\k=k(f 



(58) 



(59) 



With the help of eqs.(|58[) and ( [59| ) the integral over s of the h part of ( |5^ ) can be done obtaining 



Sex,2{R) 



n'ki r,„ -^U + 6,4 + 2fcg4^-^o4)^^'°^' ' - ' - ^' 



36m7r^ 



^{(/-l) 



(6-2feo%^^^+4^o 



/ " f 

"VfkVko 



',2' ,— + (3 + 2A;o^)AA:o 
/^ ko f 



kofk 



f f 

, (Vfco)^ 
ko 



A/ , (V/)2^ 



Afcn 



(60) 



Usung this equation and replacing the gradients of ko by gradients of p with the help of 
together with the ETF kinetic energy density given by (38), the h contribution to the exchange 
energy in the ETF approximation can be recast in the form given in eq. (|39|) . 
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FIGURE CAPTIONS 

Figure 1. Ratio of ofF-diagonal to diagonal quantal density matrix p(i?+ s/2, R — s/2)/ p[R) 
for '^^Ca as a function of the interparticle distance s for some values of the centre-of-masss 
distance R (in fm). The different curves appearing in the Figure are explained in the text. 

Figure 2. Ratio of off-diagonal to diagonal semiclassical density matrix p{R + s/2, R — 
s/2)/ p{R) for '^^Ca as a function of the interparticle distance s for some values of the centre- 
of-mass distance R (in fm). The different curves appearing in the Figure are explained in the 
text. 
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TABLE CAPTIONS 

Tablel. Extended Thomas-Fermi (ETF) Coulomb exchange energies for '^He, ^^O and '^^Ca 
compared with the quantal (QM), Slater (SL), Negele-Vautherin (NV) and Campi-Bouyssy (CB) 
values reported in ref. [0] . 

Table 2. Total quantal binding energies and root mean square radius for "^He, ^^O and '^^Ca 
obtained quantally (QM) and using the different approximations described in the text. 

Table 3. Total binding energies and root mean square radius for ^i7e, ^^O and '^^Ca obtained 
from the Strutinsky averaged method (ST) and with the different Extended Thomas-Fermi 
approaches described in the text. 
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TABLE I 





^He 


WQ 


^''Ca 


QM 


-0.86 


-2.98 


-7.46 


NV 


-0.47 


-2.31 


-6.42 


CB 


-0.78 


-2.75 


-7.03 


SL 


-0.74 


-2.75 


-7.05 


ETF 


-0.82 


-2.89 


-7.31 



TABLE II 





^He 




160 




4"Ca 






E 


< ^2 >l/2 


E 


< ^2 >l/2 


E 


< ^2 >l/2 




(MeV) 


(fm) 


(MeV) 


(fm) 


(MeV) 


(fm) 


QM 


-28.24 


1.72 


-106.62 


2.65 


-323.98 


3.36 


NV 


-23.69 


1.86 


-100.32 


2.71 


-313.49 


3.40 


CB 


-25.61 


1.80 


-101.96 


2.70 


-315.22 


3.39 


NV(KS) 


-21.75 


1.89 


-97.55 


2.73 


-309.42 


3.41 


CB(KS) 


-22.39 


1.87 


-98.27 


2.72 


-310.10 


3.40 


SL(KS) 


-17.51 


2.00 


-87.85 


2.79 


-291.78 


3.45 


ETF(KS) 


-23.78 


1.82 


-101.61 


2.69 


-317.11 


3.38 




^He 




160 




4"Ca 






E 


< ^2 >l/2 


E 


< ^2 >l/2 


E 


< ^2 >l/2 




(MeV) 


(fm) 


(McV) 


(fm) 


(McV) 


(fm) 


QM 


-29.66 


1.88 


-138.07 


2.63 


-406.47 


3.34 


NV 


-28.53 


1.91 


-135.72 


2.65 


-402.68 


3.36 


CB 


-29.04 


1.90 


-136.31 


2.64 


-403.22 


3.36 


NV(KS) 


-27.41 


1.93 


-133.65 


2.65 


-399.57 


3.36 


CB(KS) 


-27.53 


1.93 


-133.80 


2.65 


-399.62 


3.36 


SL(KS) 


-25.80 


1.96 


-129.49 


2.67 


-392.33 


3.37 


ETF(KS) 


-29.05 


1.91 


-135.38 


2.64 


-402.66 


3.35 
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TABLE III 





^He 




160 




40Ca 






E 


< ^2 >l/2 


E 


< ^2 >l/2 


E 


< ^2 >l/2 




(MeV) 


(fm) 


(MeV) 


(fm) 


(MeV) 


(fm) 


ST 


-16.71 


1.94 


-89.15 


2.75 


-296.79 


3.42 


SL 


-22.77 


1.94 


-97.73 


2.75 


-307.94 


3.42 


ETF(a) 


-22.51 


1.94 


-98.06 


2.75 


-310.51 


3.42 


ETF(b) 


-21.71 


1.94 


-96.20 


2.75 


-307.03 


3.42 


ETF{h^) 


-17.64 


1.94 


-90.19 


2.75 


-298.71 


3.42 




^He 




160 




40Ca 






E 


< ^2 >l/2 


E 


< r2 >i/2 


E 


< ^2 >l/2 




(MeV) 


(fm) 


(MeV) 


(fm) 


(MeV) 


(fm) 


ST 


-22.95 


1.95 


-125.14 


2.67 


-386.44 


3.37 


SL 


-31.85 


1.95 


-140.59 


2.67 


-409.45 


3.37 


ETF(a) 


-27.78 


1.95 


-133.39 


2.67 


-398.18 


3.37 


ETF(b) 


-27.60 


1.95 


-132.85 


2.67 


-397.22 


3.37 


ETF(?i'^) 


-23.67 


1.95 


-125.30 


2.67 


-386.07 


3.37 
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R=5 




-0.25 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



S(fm) 



R=5 




-0.25 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



S(fm) 



